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Abstract 

In the theory of Bose-condensed systems, there exists the well known problem, 
the Hohenberg-Martin dilemma of conserving versus gapless approximations. This 
dilemma is analysed and it is shown that it arises because of the internal incon- 
sistency of the standard grand ensemble, as applied to Bose-systems with broken 
global gauge symmetry. A solution of the problem is proposed, based on the notion 
of representative statistical ensembles, taking into account all constraints imposed 
on the system. A general approach for constructing representative ensembles is for- 
mulated. Applying a representative ensemble to Bose-condensed systems results in 
a completely self-consistent theory, both conserving and gapless in any approxima- 
tion. 
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1 Introduction 

Bose-Einstein condensation in Bose systems is commonly associated with the global gauge 
symmetry breaking. The most general way of the gauge symmetry breaking is by means 
of the Bogolubov shift in the field operators [1-4]. In his first works [1,2], Bogolubov 
considered a weakly-nonideal Bose gas in the limit of a small condensate depletion, which 
can occur at low temperatures and asymptotically weak interactions. The following gen- 
eralization of the Bogolubov ideas, based on the standard grand ensemble with broken 
gauge symmetry, confronted the well known problem of theories being either nonconserv- 
ing or exhibiting an unphysical gap in the spectrum. Hohenberg and Martin [5] were, 
probably, the first who analysed in detail this dilemma and classified the approaches onto 
conserving versus gapless. In a nonconserving theory, the dynamic conservation laws are 
not valid and general thermodynamic relations do not hold, which is related to the ther- 
modynamic instability of the system. If the theory is conserving, but exhibits a gap in 
the spectrum, it also has incorrect thermodynamics and contradicts the theorems of the 
gapless spectrum by Bogolubov [4] and Hugenholtz and Pines [6]. A thorough discussion 
of the Hohenberg-Martin dilemma has recently been done by Andersen [7], where the 
extensive literature on the problem is cited. 
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It is important to emphasize that the inconsistency, happening in the theory, is in- 
ternal, but not connected with the employed approximations. Really, suppose one ap- 
plies the standard grand ensemble, with a chemical potential fi, to an equilibrium Bose- 
condensed system. After breaking the gauge symmetry, by means of the Bogolubov shift, 
one has to minimize the grand thermodynamic potential with respect to the number of 
condensed particles A^o- Bogolubov [1~4] illustrated this procedure for a weakly nonideal 
gas. Ginibrc [8] proved that this procedure is general for arbitrary Bose systems. The 
Bogolubov- Ginibrc minimization procedure defines the chemical potential /j, I^Boip, T) 
as a function of density p and temperature T. 

On the other hand, there exists the Hugenholtz- Pines relation [6] also defining the 
chemical potential p — > ij,hp{p,T). The most general proof on this relation, for any 
Bose system and arbitrary temperatures, is based on the consideration of increments of 
statistical averages under gauge variations [4]. This is analogous to the consideration of 
the Ward-Takahashi identity. The Hugenholtz-Pines relation follows from the symmetry 
of the Hamiltonian, but is not connected with the system stability. 

Thus, there are two representations for the chemical potential. One, /ibg, is a conse- 
quence of the condition of thermodynamic stability, while another, fiHP, results from the 
symmetry properties. Two different, both physically and mathematically, conditions do 
not need to define one and the same quantity, so that generally, Pbg ^ Php- And this is 
exactly what happens [5,7], making the theory not self-consistent. 

So, it is not some particular approximations that are guilty for the arising inconsis- 
tency. Also, the gauge symmetry breaking is not to be blamed. The concept of symmetry 
breaking under phase transitions is rather general [9,10]. But the inconsistency, related 
to the inequality Hbg Php, is generic for the used statistical ensemble. 

The Hohenberg-Martin dilemma has been addressed in a number of papers. One 
usually attempts to cure the problem in a mean-field approximation complemented by 
some phenomenological tricks. The most common, and the least justified, trick is by 
omitting the anomalous averages. One often ascribes this trick to Popov, terming it 
the "Popov approximation". First of all, however, it is not difficult to infer from the 
original works by Popov [11-13], cited in this regard, that he has never suggested or used 
such an unjustified trick. The person who really first proposed and employed the trick 
was Shohno [14]. This was known in literature as Shohno model. But in the case of 
the broken gauge symmetry the anomalous averages become principally important and 
cannot be omitted. Direct calculations [15] show that the anomalous averages at low 
temperatures can be even much larger than the normal ones. So, neglecting the terms 
that are larger than what is left cannot be a reasonable mathematical approximation. 
In addition, the omission of anomalous averages renders the system thermodynamically 
unstable [15,16], with incorrect thermodynamics [17]. 

To remove the gap in the spectrum, one sometimes resorts to what Bogolubov [4] 
called the mismatch of approximations. Then one defines the chemical potential /ibg 
from the condition of thermodynamic stability in one approximation, but, for cancelling 
the appearing gap, one invokes the self-energies in another, higher, approximation, or 
just adding phenomenological terms. As is evident, the mismatch of approximations is 
not a logical procedure, is not uniquely defined, and does not solve the problem that 
for each given approximation the theory remains not self-consistent. Also, the mismatch 
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of approximations leads to incorrect thermodynamics, with a first-order phase transition, 
instead of the second-order one [18-20]. Thermodynamic self-consistency can be preserved 
only when all quantities are derived from the same Hamiltonian in the same chosen 
approximation [21,22]. 

As is mentioned above, the inconsistency, manifesting itself in the inequality Hbg 
IJ'HP, is generic for the used statistical ensemble and does not depend on the employed 
approximations. The problem can be solved by using a representative statistical ensemble, 
which takes into account all constraints imposed on the Bosc-condensed system [23] . Then, 
for instance, the Hartree-Fock-Bogolubov (HFB) approximation can be made conserving 
and gapless [24]. 

In the present paper, the general notion of representative statistical ensembles is formu- 
lated for arbitrary statistical systems. The representative ensemble for a Bose-condensed 
system is defined. It is shown that the theory, based on the representative ensemble, 
is completely self-consistent, independently of the used approximations. This solves the 
Hohenberg- Martin dilemma, suggesting an approach that is free from inconsistencies and 
paradoxes. 

Everywhere, below the system of units is used, where h = 1 and ks = 1- 

2 Representative Ensembles 

It was already Gibbs [25] himself who mentioned that for a correct description of a sta- 
tistical system, in addition to a canonical distribution, one has to accurately take into 
account all conditions and constraints, imposed on the system and uniquely representing 
its properties. From this general idea stems the term "representative ensembles", which 
was used by ter Haar [26,27]. For equilibrium systems, one defines representative ensem- 
bles through the conditional maximization of the Gibbs entropy, as was done by Janes 
[28]. A generalization for quasi-equilibrium systems was given in Ref. [29]. Here, we 
describe a general method for constructing representative ensembles for both equilibrium 
and nonequilibrium systems. 

An equilibrium statistical ensemble is a pair {JF, p} of a space of microstates JF and a 
statistical operator p on this space. Let a set {Q} of self-adjoint operators Q be given, 
which we shall call the condition operators, since these define the statistical conditions 

C, = <C,>= Trp A , (1) 

which are assumed to be imposed on the considered system. The trace operation in 
Eq. (1) is over the space of microstates J-. One obvious constraint is the normalization 
condition for the statistical operator, 

1 = < > = Trp , (2) 

where Ijf is the unity operator in !F. Another standard condition is the definition of the 
internal energy 

E = < H > = Tip H (3) 

as the average of a Hamiltonian H. However, in addition to these standard conditions, 
other statistical conditions (1) can also be necessary for correctly representing the system. 
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For the conditional maximization of the Gibbs entropy S = — Trpln p, under statistical 
conditions (1), (2), and (3), it is convenient to introduce [30] the information functional 

7[p] = IVp In p + Ao (Trp - 1) + /3 (Xrp^ - e) + P^.^'i {^^P^i - Ci) , (4) 

i 

with the Lagrange multipliers Aq, P, and The statistical operator is given by the 
minimum of the information functional (4) , 

P = \ e-^^ , (5) 

with the grand Hamiltonian 

H = H + J2uA. (6) 

i 

Here Z = exp(l + Aq) is the partition funtion and (3 = T"^ is inverse temperature. 

The representative statistical ensemble for an equilibrium system under constraints (1), 
(2), and (3), is the pair {J^, p} of a space of microstates and the statistical operator 
(5) with the grand Hamiltonian (6). 

Note that the condition operators Ci are to be self-adjoint, but they do not need to 
be compulsorily the integrals of motion. Thus, the number-of-particle operator TV is not 
an integral of motion, when the gauge symmetry is broken. 

A nonequilibrium statistical ensemble is a triplet {JF, p, d/dt} of a space of microstates 

the statistical operator p = p(0) at the initial time t — 0, and of a prescribed evolution 
law, symbohcally denoted here as d/dt. The law of evolution can be given as the Liouville 
equation for the time- dependent statistical operator p(t) or as the Heisenberg equations 
of motion for physical operators, or in some other form. A general way of deriving the 
evolution equations is by extremizing an effective action functional with respect to the 
given dynamical variables [31]. 

To construct an action functional, one defines the corresponding Lagrangian 

L = E-H , (7) 

in which E is the time-dependent energy operator, containing the time derivative. The 
effective action, for a system with the Lagrangian (7), under constraints (1), is 

A = j (L- uA) dt , (8) 

where Vi are the Lagrange multipliers. Equation (8) can be rewritten as 

A= I (E-H)dt, (9) 



with the same grand Hamiltonian (6). The evolution equations follow from the extrem- 
ization of the effective action, 6 A = 0. For instance, if the effective action A = A[-(/'] is a 
functional of field operators ip{r,t) and ip'^{r,t), then the latter satisfy the equations 

sm mL^o. (10) 



Sij^{r,t) ' Sijj{r,t) 

The grand Hamiltonian (6) governs both the statistical properties of equilibrium systems 
and the temporal evolution of nonequilibrium systems. 
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3 Bose Systems 



Let us consider a one-component system of spinless particles characterized by the field 
operators t/j and satisfying the Bose commutation relations. For the system above the 
critical temperature Tc, the space of microstates is the Fock space ^(^) generated by 
the field operators [30-32]. In addition to statistical conditions (2) and (3), there is one 
normalization condition for the total number of particles N =< N >. Hence, the grand 
Hamiltonian (6) is H = H — /iN, where /i is the usual chemical potential. 

However, the situation changes below the critical temperature T^., when the global 
gauge symmetry becomes broken. The symmetry breaking is realized by the Bogolubov 
shift 

^{r,t) 4^{r,t) = 77(r,i) + Vi(r,t) , (11) 

in which ri{r, t) is the condensate wave function and i/ji{r, t) is the field operator of uncon- 
densed particles. The field operators ipi and ipl generate the Fock space J-'{tpi), which is 
the space of microstates for the system with broken gauge symmetry. The spaces of mi- 
crostates ^(^), for the system above and for the system below are mutually 
orthogonal [33,34]. 

For the system with broken gauge symmetry, there are now two variables, the conden- 
sate wave function rj and the field operator of uncondensed particles ipi. These variables 
are assumed to be linearly independent, being orthogonal to each other, 

Jrj*{r,t)M^,t)dr^O. (12) 

Respectively, instead of one normalization condition for the total number of particles, we 
have now two normalization conditions, for the number of condensed particles 

TVo = < A^o > , (13) 

where, the average is over the space J-'{tjji), Nq = Noljr^ ijr being a unity operator in 
J^{i>i), and 

No^J \v{r,t)\'dr, (14) 
and for the number of uncondensed particles 

Ni = < TVi > , (15) 

with the operator of uncondensed particles 

N, = J jjjl(r,t)M-^,t) dr . (16) 

The total number of particles is = A^^q + A^i . 

With the broken gauge symmetry, the average < '^i > might be nonzero. This, 
however, would mean that quantum numbers, such as spin or momentum, would not be 
conserved. To avoid this, one implies the constraint 

<V^i(r,t) >= 0, (17) 
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guaranteeing the conservation of quantum numbers. Constraint (17) can be represented 
in the standard form (1) of a statistical condition by introducing the self-adjoint operator 



A = y [A(r, t)4{r, t) + A*(r, t)^i(r, t)] dv , (18) 



in which A(r, i) is a complex function. Then constraint (17) is equivalent to the quantum 
conservation condition 

< A > = . (19) 

Thus, for a system with broken gauge symmetry, there are three statistical conditions, 
(13), (15) and (19). To construct the corresponding representative ensemble, we have to 
follow the general prescription of Sec. 2. 

For a system in equilibrium, in view of Eq. (4), the information functional is 

I[p] = Trp In p + Ao (Trp - 1) + /9 (TxpH - E) - 

-/3/xo (TrpiVo - A^o) - /i/xi (TrpiVi - A^i) - /3Tr/)A . (20) 

The minimization of this functional yields the statistical operator (5) with the grand 
Hamiltonian 

H^H- pioNo - f^iNi - A , (21) 

in agreement with Eq. (6). 

In the general case of a nonequilibrium system, we need to consider the Lagrangian 
(7), in which the temporal energy operator is 



E 



I 



rj*{r, t) i ^ 7;(r, t) + 4{r, t) ^ ^ V'l (r, t) 



dv . 



The action functional (9) has to be extremized with respect to two linearly independent 
variables, 77 and ■^i, which implies 

67j*{r,t) ' S4{r,t) ■ ^ ^ 

The latter is equivalent to the equations 

■ 9 , SH[r],ipi\ ■ 9 , . SH[r],ipi\ 

I — mv, t) = — , I — ipiir, t) = — I . (23) 

dt ' ' 5r]*{r,t) dt ^ SiPl{r,t) ^ ' 

The average of the grand Hamiltonian (21) is 

<H>^ E- hqNq - iiiNi . (24) 

If in experiment only the total number of particles, or the total average density, can be 
fixed, then average (24) should be expressed in the standard form < H >— E — jiN 
through the system chemical potential Comparing this with Eq. (24) defines the 
chemical potential 

H = HqUq + jjiiUi (25) 
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through the Lagrange multiphers /lo and /ii and through the fractions uq = Nq/N and 
m = Ni/N. 

Let us take the Hamiltonian energy operator in the usual form 

H = j ^t(r) ^Y_ + u^ i>{Y) dr + 

+ ^ J V'^(r)V^^(r')$(r - r')V'(r')V'(r) drdr' , (26) 

in which ■0(r) = il^{r,t) is the Bogolubov shifted field operator i/j = rj + t/ji, defined in 
Eq. (11). Here U = U{r,t) is an external potential. The interaction potential $(r) is 
symmetric, so that $(— r) = $(r), and is assumed to possess the Fourier transform. With 
the Bogolubov shift (11), Hamiltonian (26) can be represented as a sum of five terms, 
H = + + ^ ^(3) ^ ^(4)^ depending on the powers of the operators t/ji and 
ipi in each of them [34]. The zero-order term contains no ipi. The necessary and 
sufficient condition for the validity of constraint (19) is the absence in the Hamiltonian of 
terms linear in -01 [35]. This requires to define the Lagrange multipliers A(r, t) such that 
H('^ = 0. 

The evolution equations (23) give 

i ^ ri{v, t)^(^-^ + U-^io^ 7^(r, t) + J $(r - r') [|,7(r') |^77(r) + X(v, v')] dr , 



+ 1 $(r - r') [|r7(r')| Vi(r) + r/*(r')r/(r)V'i(r') + r/(r')r/(r)^I(r') + X(r, r')] dr' , (27) 

where, for brevity, the time dependence in the integrands is not shown explicitly and the 
notation for the correlation operator 

X{r,r')=^l{r')Mr'Hr)+ 

+^I(r')r/(r')^i(r) +r/*(r')t^i(r')^i(r) +^I(r')^i(r')^i(r) 

is used. It is straightforward to check that the evolution equations (27) guarantee the 
validity of all local conservation laws on the operator level, as it should be for the equations 
derived from a variational procedure. 

Averaging the first of Eqs. (27) over ^{ipi), we obtain the explicit equation for the 
condensate wave function 

d / \ 

+ 1 $(r - r') [p{r')r]{r) + pi(r, r')r/(r') + ai(r, r')r]*{r') + ^(r, r')] dr' , (28) 

in which 

pi(r,r') = < ipl{r')i!i{r) > , ai{r,r') = < iPi{r')iPi{r) > , 



7 



p(r) = po(r) +pi(r) , po(r) = |?7(r)|' , Pi(r) = pi(r,r) = < V'i(r)V'i(r) > . 

Equation (28) is valid for an arbitrary Bose system. For a system in equilibrium, we 
require that the condensate function be time independent, drj/dt — 0. If, in addition, the 
system is uniform, with no external fields, then Eq. (28) yields 



//o = p$o + / *(r) 



Pi(r,0)+<7i(r,0) + ^^'"'°^ 



dr . (29) 



The matrix Green function ^(12) for uncondensed particles is defined in the usual 
way [4-6], where the common abbreviation can be employed, denoting the set {rj,tj} by 
a single number j. The elements of this matrix are 

G'n(12) = -i < ^^^1(1)^1(2) > , 6-12(12) = -i < tV'i(l)V'i(2) > , 

621(12) = -i < f4{l)4{2) > , 622(12) = -i < f4{l)M'2) > , (30) 

where T is the chronological operator. For an equilibrium system, we may pass to the 
Fourier transform G(k, a;) of G(12). Equations of motion give 

a; + A;V2m + Sii(k,a;) -/xi 

Gi2(k,^) = - = G2i(k,a;) , (31) 

where E(k, a;) = [Ecj^(k, a;)] is the self-energy matrix, and the denominator is given by 
the equations 

D{k,uj) ^ S2,(k,u;) -Gr/(k,a;)Gr/(k, -a;) , 

G^/(k,a;)=a; - ^ - Eii(k, u;) + //i . 
From the Bogolubov theorem [4], we have 

|Gii(k,0)-Gi2(k,0)|>^. (32) 

Using here Eqs. (31), we obtain the Hugenholtz- Pines relation 

//i = Eii(0,0)-Ei2(0,0) , (33) 

which makes the spectrun gapless. 

Comparing Eqs. (29) and (33), we see that these expressions are rather different, and 
there is no reason to require that /Iq be equal to Hi. As an illustration, we can resort to 
the Hartree-Fock-Bogolubov (HFB) approximation, in which we find 

1 1 

^ p^o ^ p^O 
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where V is the system volume, $p is the Fourier transform of $(r) and Up and ap are the 
normal and anomalous averages, respectively, in the momentum space. 



with ttk being the Fourier transform of ipi{r). Clearly, ^ ^i. The anomalous average Ck, 
as can be shown by direct calculations [15,24], can be much larger than the normal average 
Uk- Hence, there is no any sense in omitting a^. As has also been mentioned above, the 
presence of cr^ is crucially important for the stability of the system [15,16,23,36]. 

The physics of dilute Bose gases has resently attracted a great deal of attention, both 
experimentally and theoretically (see review works [7,37-40]). These gases are described 
by the local interaction potential 



$(r) = %5{r) , 



$0 = 47r 



m 



in which is the scattering length. For this potential, Eqs. (34) reduce to 
A^o = (P + Pi + o-i)'^'o , /ii = (p + pi - ai)$o , 

where 



(35) 



(36) 



1 



Pi = y Y.'^k, 



V 

^ k^Q 



The particle spectrum has the Bogolubov form 



£k 



\ 



with the sound velocity c defined by the equation 

mc^ = (po + o-i)$o , 



(37) 



where po = Nq/V is the condensate density. The anomalous average ai, invoking the 
dimensional regularization, can be represented as 



(mc)^ / — (mcf /■= 

TT^ ^ 2\/27r ^0 

The condensate fraction becomes 



(Vl 



1)1/2 



Vi 



coth 



mc 
~2T 



- 1 



dx . (38) 



no 



1 



(mc)^ 
37r^p 



3 / \ 



1/2 



coth 



mc 



X 



dx , 



(39) 



And for the superfiuid fraction, defined through the response to a velocity boost [37], we 
obtain 

[mcj^mc^ (Vl + x"^ — l)^^'^x dx 



eV^T^^pT Jo y/l + x^smh'^(mc'^x/2T) 
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The above equations show that the system undergoes the second-order phase transition 
at the critical temperature 



T 

c 



271 

m 



P 



2/3 



.C(3/2). 

coinciding with that for the ideal gas, as it should be for a mean-field approximation with 
the local interaction potential (35). Note that all integrals in Eqs. (38), (39), and (40) 
are well defined and convergent. 

In conclusion, the old-standing problem, the Hohenberg- Martin dilemma, is analysed, 
which classifies the theories of Bose-condensed systems onto either conserving or gapless. 
This dilemma is shown to be of general nature, not dependent on the used approximations, 
and having its footing in the generic inconsistency of the standard grand ensemble applied 
to the Bose systems with broken global gauge symmetry. The notion of representative sta- 
tistical ensembles for arbitrary systems is formulated. A representative ensemble uniquely 
represents the given statistical system by taking into account all constraints and condi- 
tions imposed on the system. The usage of representative ensembles for Bose-condensed 
system provides the remedy from internal inconsistencies, cures old paradoxes, and gives 
the solution to the Hohenberg-Martin dilemma. The suggested theory for Bose-condensed 
systems, based on a representative ensemble, is completely self-consistent, conserving, and 
gapless. This is achieved by introducing two Lagrange multipliers, /xq and /xi, which are 
responsible for two normalization conditions, for Nq and A^i, respectively. Note that in the 
current experiments with trapped atoms [37-40] two atomic numbers can be controlled, 
the total number of atoms N and the number of condensed atoms A^o- In addition, the 
operator Nq is the integral of motion, since [TVq, H] — 0. So, fixing two normalization 
conditions is justified both mathematically as well as physically. Two Lagrange multi- 
pliers are necessary in order to uniquely define the condensate wave function, which is 
a two-component order parameter. This agrees with the general rule for any statistical 
system: 

To uniquely define an order parameter of a system, it is necessary to introduce as 
many Lagrange multipliers as is the number of components in the order parameter. 
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